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Based on the recently proposed band model, the electronic specific heat of moderately heavy 
electron compound YbAl3 are investigated. The band term of the Hamiltonian consists of three 
parts; conduction electrons described by the nearly free electron method, localized 4f electrons 
of Yb ions and the hybridization term between these electrons. Extracting several bands near 
the Fermi level, we reconstruct the low-energy effective Hamiltonian in order to consider the 
correlation effect, which is studied by using the self-consistent second order perturbation theory 
combined with local approximation. The temperature dependence of the specific heat c v (T) is 
calculated as a function of temperature T from the numerical derivative of the internal energy. 
Sommerfeld coefficient 7 is also calculated from the direct formula. The overall structure of 
Cv(T)/T is in quantitative agreement with the experimental results, which have the character- 
istic two-peak structures. They originate from the correlation effect and the structure of the 
non-interacting density of states, respectively. We show that our effective Hamiltonian yield- 
ing the realistic band structure may describe quantitatively heavy electron compounds with 
conduction bands composed of s- or p- electrons. 

KEYWORDS: YbAl3, nearly free electron model, empty-core pseudo-potential 
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1. Introduction 

The periodic Anderson model (PAM) and its exten- 
sions are often employed in order to investigate the physi- 
cal properties of heavy fermion compounds. 1 These mod- 
els succeeded in capturing the qualitative properties. 
However, much simplification is done in PAM, sometimes 
leading to incorrect results even qualitatively, and it is 
difficult to describe quantitative properties of realistic 
materials. Therefore, tight-binding model is constructed 
form the band calculations based on the local density 
approximation (LDA) and correlation effects are incor- 
porated in various ways. 

Recently we proposed a new method of the construc- 
tion of the effective Hamiltonian with the realistic band 
structure for heavy fermion materials whose conduction 
bands consist of s- and p- electrons mainly. 2 By using 
the combination of solid methods, we tried to construct 
the quantitatively reliable model with only a few pa- 
rameters. It is well known that the band calculation for 
heavy fermion compounds sometimes reproduces incor- 
rect f electron level. Thus, we can study the realistic 
materials from another point of view by means of this 
technique. Furthermore, this model can capture the ba- 
sic properties of the realistic materials such as the lat- 
tice symmetry, and the various physical quantities can 
be calculated quantitatively from almost absolute zero 
to finite temperature region within the desired accuracy 
in the wave- vector k and frequency u> spaces. Thus, this 
model may be good starting point to discuss the physical 
properties of strongly correlated systems. 

The Hamiltonian consists of conduction electrons de- 
scribed by the nearly free electron (NFE) method, lo- 
calized f electrons and the hybridization between the 
conduction and the localized f electrons. This method 
includes only a few parameters. We apply it to typi- 



cal heavy fermion compound YbAl3, which has NFE- 
like conduction bands below and around the Fermi en- 
ergy Ep. The obtained band structure is in good agree- 
ment with the result of the band calculation provided 
by Harima, 2,3 and the optical conductivity is also con- 
sistent with the experimental result. 4 However, the cal- 
culated value of Sommerfeld coefficient 7 is about 13.5 
mJ-K~ 2 -mol -1 , 5,6 which is 1/3 of the experimental result 
40 mJ-K _2 -mol -1 . 6,7 This discrepancy may result from 
the lack of the correlation effect which was not taken into 
account in the previous study. 

Note here the physical properties for Yb AI3 . The mean 
valence is in the range from 2.65 to 2.8 s ' 9 and the Som- 
merfeld coefficient 7 is about 40 mJ- K -2 - mol -1 , which 
indicates this material is a valence fluctuation compound 
and does not have strong correlation effect. The ampli- 
tude of the crystal field splitting of 4f J = 7/2 level is 
below 10 meV, 10 which is much smaller than the Kondo 
temperature Xk ~ 500 K. 6 In the low-temperature re- 
gion, the specific heat divided by temperature c v (T)/T 
shows prominent temperature dependence; a peak struc- 
ture appears at T ~ 70 K where the amplitude is about 
80 mJ- K~ 2 - mol -1 . Furthermore, c v (T)/T decreases 
from T = below T ~ 25 K, which was interpreted as 
the effect of the lattice coherence. 6 The coherence tem- 
perature is about Tech ~ 35 K, below which the electrical 
resistivity is proportional to T 2 . 6,11 This result indicates 
that YbAl3 satisfies the Fermi liquid behavior. The op- 
tical conductivity shows the mid-infrared peak structure 
at hw ~ 0.25 eV. 4 ' 12 

In this paper, we focus on the heavy fermion metal 
YbAl3 and construct the low-energy effective Hamilto- 
nian with the realistic band structure and the Coulomb 
interaction. The latter term is not included in the pre- 
viously proposed effective Hamiltonian. We investigate 
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The plane wave with wave vector k and the localized 
4f-electron can be expanded by the spherical harmonics, 
so that the matrix element of the hybridization term is 
also described by this spherical harmonics. For simplicity, 
we neglect k-dependence of the radial part and regard 
this part as a constant which denotes the hybridization 
amplitude V. We stress that although we neglect the 
k-dependence of the radial part, we keep the angular 
dependence of k completely. 

The constructed band term of the Hamiltonian is given 
by 



^band 



/ # T C o 

HI 



H hyb 



(3) 



Diagonalizing the band term of the Hamiltonian, we can 
rewrite as 

H hand = J2 E ^ai a a ka , (4) 



where E^ a denotes an energy eigenvalue, and a^a is an 
annihilation operator of a quasi-particle with the wave- 
vector k and the band index a. The energy bands are 
shown in Fig.l, which is in rough agreement with the 
LDA result around and below the Fermi level except M- 
point. 3 The topology of the Fermi surface is also consis- 
tent with the LDA result except the same point. 

Next, let us consider the interaction term H mt . The 
correlation effect may be essential on the localized 4f- 
electrons at the low-temperature region, so we assume 
that H lnt is given by 

H int = U Yl riUMnuM', (5) 

i,M<M',7=| 

where riiju is a number operator of the 4f electron with 
J. M is the z-component of J and U stands for the 
Coulomb interaction. 

By using the diagonalizcd bases for the band term, 
an annihilation operator of the 4f electron with the 
wave-vector k and the angular momentum J and its z- 
component M is written as 



the physical properties and discuss the efficiency of the 
Hamiltonian. Applying the self-consistent second order 
perturbation theory with local approximation to this 
low-energy effective Hamiltonian, we calculate c v (T)/T 
and the Sommerfeld coefficient. 

This paper is organized as follow. We review the band 
model for YbAl3 with the realistic band structure and 
construct the low-energy effective Hamiltonian with the 
Coulomb interaction to discuss the correlation effect in 
the next section. In §3, we show the density of states 
(DOS), the Sommerfeld coefficient and the specific heat 
divided by temperature. Summary and discussions are 
given by in §4. 

2. Effective Hamiltonian for YbAl 3 and Method 

In this section, the low-energy effective Hamiltonian 
for YbAl3 is constructed, which can be written as follows 

H cff = ^band +H int^ ^ 

iJ band = H C + H f + iJ hyb , (2) 

where H hand (H lnt ) represents the band (interaction) 
part, respectively. H hand consists of the conduction band 
H c , the localized 4f orbitals H l , and the hybridization 
between the conduction and the localized 4f-electrons 
term _ff hyb , while H mt consists of the Coulomb inter- 
action between the localized 4f-electrons of Yb ions. 

We briefly review the band part of the Hamilto- 
nian H hand based on the proposed method. 2 Each elec- 
tronic configuration of YbAl 3 is [Ne](3s) 2 (3p) 1 for Al and 
[Xc](4f) 14 (5d)°(6s) 2 for Yb atom. There are 25 electrons 
per unit cell of YbAl 3 except the closed shells. The mean 
configuration of the localized 4f-electrons lies between f 14 
and f 13 , which indicates that YbAl 3 shows the valence 
fluctuation behavior and becomes a typical metallic sys- 
tem. Most of the conduction electrons consist of s- and 
p-electrons, which wholly spread over the crystal with the 
weak periodic potentials. Therefore, for the conduction 
bands, we employ the NFE method with the empty-core 
pseudo-potential, whose screening effect is considered by 
using the dielectric function of the free electron gas. The 
atomic potential and the core radius in each atom are 
employed from Harrison's textbook. 13, 14 

The 14-fold degenerated 4f-electron levels split into the 
total angular momentum J=7/2 (8-fold degeneracy) and 
J=5/2 (6- fold degeneracy) levels because of the strong 
spin-orbit coupling in the case of YbAl 3 . The difference 
between both energy levels is of the order of 0.1 Ry, which 
is obtained from the band calculation. 3 The J=7/2 level 
Ej^ 2 is located near and below the Fermi level, while 
the J=5/2 level is far from the Fermi level. 5 Although 
the J=5/2 level may not affect the physical properties in 
the low-temperature region, we also include J=5/2 term 
in the Hamiltonian iJ band and set the energy difference 

7 /2 

between both levels as 0.1 Ry. Thus E { is a tuning 
parameter in this method. Note here that the amplitude 
of the crystal field splitting of 4f J = 7/2 level is very 
small compared to the Kondo temperature. Therefore, 
although to consider the crystal field splitting effect is 
easy in this framework, it is not discussed in this paper. 



/kJM = ^WkJMQflka, (6) 

a 

where UkjMa stands for the eigenvector of JJ band . Thus, 
H lnt is written as 

" int = ^E E U^ai + ^ a ,at,_^av ,, (7) 

k.k'.q aa'/3/3' 

where N is the number of the lattice sites, and the inter- 
action matrix cT^ k'^q 1S as f°ll° ws 

^k.k'fq = U E W k+qJMa<JMa' u k'-qJM'/3 u k'JM'/3'- (8) 
M<M' 

The hybridization does not affect the band disper- 
sions drastically. There are almost fiat bands near and 
below the Fermi level, which indicates that these diago- 
nalizcd bands still hold the 4f character strongly (shown 
in Fig.l). Thus we extract 14 bands near the Fermi level, 
which almost lie in \E\^ a — Ep\ < 0.1 Ry and include 
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Fig. 1. Energy band dispersion B kce of YbAb. The 4f J = 7/2 

7/2 

level is set at B f = 0.6 Ry. The solid lines stand for the ex- 
tracted bands in order to consider the correlation effect. The 
dashed lines are the bands not included in the effective Hamil- 
tonian eq. (10). The dot-dashed line stands for the Fermi level. 
The lower panel is the magnification of the upper panel near the 
Fermi level. 



all the 4f J = 7/2 components. There are sufficient to 
discuss the effect of the Coulomb interaction. 

To avoid the extremely complicated treatment, we re- 
gard the interaction matrix as diagonal within the ex- 
tracted band indices as in eq. (9) below. Then, the inter- 
action term can be expressed in the Hubbard type with 
the number operators of the extracted bands, n^n^. In 
addition, the amplitude is approximated by a constant. 
Thus the interaction matrix is written as 



U8 nn '8r, 



(9) 



This treatment means that we allow Coulomb interaction 
to act both on the 4f and conduction electrons equally. 
However, it acts on the 4f electrons more strongly since 
the DOS is much larger. Hereafter, we regard the ef- 
fective Coulomb interaction U as a parameter to dis- 
cuss the correlation effect. Since Coulomb interaction 
is strongly screened and renormalized in metallic heavy 
fermion compounds, the amplitude may not be so large 
at the low-energy region. 

Note that the band index is labeled in the order of 
energy eigenvalue from smallest to largest at each wave- 
vector k, for simplicity. One should classify the band in- 



dex according to the space symmetry, but the simplifica- 
tion for the labeling is irrelevant to the physical quanti- 
ties such as specific heat in the present study. 

Finally, we obtain the following low-energy effective 
Hamiltonian 



H 



off 



= E^ 

ka 



ka a ]La a ka 



+ 



-Y, 

k,k',qa</3 



k +qa aka4'-q/? a k'/3- (10) 



Here, Ek a consists of only the extracted bands near the 
Fermi level. 

Let us consider the correlation effect based on this 
low-energy effective Hamiltonian with the realistic band 
structure. The amplitude of the Coulomb interaction is 
not so large because of the renormalization effect at the 
low-energy region, so that we employ the self-consistent 
second order perturbation theory with respect to the 
Coulomb interaction. Furthermore, since the 4f electrons 
have strongly localized character, local approximation is 
applied. Note that the off-diagonal components in the 
self energy vanish in the present approximation since the 
Hubbard-type interaction is employed. 

The Green function and the local Green function are 
defined as 

G a (k,iw„) = [iu>„ + fi - E ka - E a (iw„)]~ ,(11) 

gL oc m = ^E G «( k >«)> ( 12 ) 

k 

where u n stands for Matsubara frequency and \i is the 
chemical potential. Combined with local approximation, 
the self energy is given by 

E a (iw„) = UT J2 G£ c (W)e iw "'° + 



-U 2 T 2 



E 



xG^(icj m + \w m ,)G^(w m ,), (13) 

where + is a positive infinitesimal number. T denotes 
temperature which includes the Boltzmann constant and 
has a dimension of energy. The internal energy is written 
as 



E" 



ka I n ) 

NT 

+ -^E G « c ( iw ") s «( iw «)- 



(14) 



The sum of Matsubara frequency is replaced by the in- 
tegration of real frequency by means of analytic contin- 
uation. 

The specific heat is obtained from the derivative of 
the internal energy with respect to temperature, which 
is given by 

iV A d£ int 



Cv(T) 



N AT ' 



(15) 
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where N\ denotes Avogadro constant. The estimation 
of the derivative is replaced actually with the numerical 
difference. 

The Sommerfcld coefficient 7 is given by 



3000 



7 



3N 



kci' 



1 



OJ = 



A a (k, 



-Im G a (k,uj), 



16) 



(17) 



within the Fermi liquid theory with local approximation. 
j4 Q (k, 0) stands for the spectral function on the Fermi en- 
ergy. For metallic systems with the Fermi liquid behavior, 
limj^o c v(T)/T must be equal to 7. Note that YbAl 3 ex- 
hibits the Fermi liquid property at the low temperature 
region (T < 35 K). 

The total DOS and the number of 4f-electrons with 
the angular momentum J are given by 



(18) 



uwjMa\ 2 J dcu A a (k, w)/(w%19) 



1 

TV 



E 

k,M.Q 



where f(to) stands for Fermi distribution function. 
In the present method, the Hamiltonian includes only 

7 /2 

a few parameters (4f J = 7/2 level E f , the hybridiza- 
tion amplitude V and the effective Coulomb interaction 
U). Compared with some experimental results, these pa- 
rameters are estimated. 

3. Results 

In the NFE method, the weak periodic potential is 
given by the sum of reciprocal lattice vectors as 

3 iG r r 



V(T) = £V Gl 



Gi — ^Gq + / 2 Gq + / 3 Gq. 



(20) 



(21) 

Gg (/j, — 1, 2, 3) stands for the primitive translation vec- 
tor of the reciprocal lattice and I denotes an integer set 
(= (li, I2, 13)) within our target energy range. The calcu- 
lation, however, converges within |Z^| < 3 (fj, = 1,2,3). 
We divide the first Brillouin zone into 80 x 80 x 80 pieces 
to carry out wave vector k sums. In order to focus on 
the correlation effect, the tuning parameters for the band 
structure are given by the same as those in the previous 
paper. 2 Namely, the 4f J = 7/2 level and the hybridiza- 
tion amplitude are fixed mainly as E 1 ^ 2 = 0.6 Ry and 
V = 0.03 Ry, respectively. The 4f J = 5/2 level is defined 



as E, 



5/2 
f 



E 



7/2 
{ 



0.1 Ry. The pseudo-potential parame- 
ters are defined as follows; core radius rj b = 2.0 A for 
r ' = 0.61 A for Al. 13 ' 14 Note that a larger 



Yb and - A1 



value for rj h than that of ref. 13 is chosen. 2 The valence 
of each atom is Z Yh = +2 and Z A1 = +3, respectively. 

Figure 2 shows the total DOS for a few choices of U 
near the low-energy region and at the low-temperature 
region. The prominent peak structures appear below the 
Fermi level (u = 0), which mainly consist of the 4f J — 
7/2 electrons. Although the spectrum is broadened due 
to the hybridization effect, the localized character of 4f 
electrons still remains. 
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Fig. 2. DOS as a function of the effective Coulomb interaction 
U for E 7 { /2 = 0.6 Ry at T = 2.0 X 10~ 5 Ry. The Fermi level is 
set to u> = 0. The inset shows the magnification near the Fermi 
level. 



The simple PAM or its extension produces a hybridiza- 
tion gap easily even in the metallic system, which is 
caused by the simplification, such as neglect of the wave 
vector dependence in the hybridization term. In the 
present model, the hybridization still keeps the angular 
dependence of k, so that a hybridization gap never ap- 
pears. Figure 2 shows only a pseudo gap structure. Note 
that the Coulomb interaction does not affect topology of 
the Fermi surface within the present parameter region. 

With an increase of the effective Coulomb interaction, 
the peak structure is broadened and is shifted to the 
Fermi level gradually, which means that the effective 4f 
J = 7/2 level is shifted by U while the amplitude at 
u) ~ seems almost insensitive (Fig. 2 inset). Although 
the number of 4f electrons decreases due to this shift, 
the deviation from U = Ry is not serious, which is 
shown in Fig. 3 for Ej' 2 = 0.6 Ry (and other parameters 
for reference). The difference between the case of U = 



13.6 




0.01 
U(Ry) 



Fig. 3. The number of 4f electrons as a function of the effective 
- 7/2 
Coulomb interaction U for a few choices of E { at T = 2.0 X 

10- 5 Ry. 
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Fig. 4. The Sommerfcld coefficient as a function of the effective 
Coulomb interaction U for a few choices of Ej^ 2 at T = 2.0 X 
fCT 5 Ry. 



Fig. 5. c v (T) /T as a function of temperature T for several choices 
of U. The inset shows the temperature dependence of c v (T). The 
arrows stand for the data of the Sommerfeld coefficients at T ~ 
where each symbol corresponds to that of c v (T)/T. 



Ry and that of U ~ 0.0175 Ry is less than 0.1, which 
still stays in the following experimental result. The mean 
valence of Yb ion in YbAl 3 is 2.65 ~ 2.8, which indicates 
the number of 4f electrons is the range from 13.2 to 13.35. 
Of course, further increase of U pushes up the effective 
4f level above Ep, then m decreases quickly. 

The effective Coulomb interaction is introduced for the 
diagonalizcd bands which partially include the conduc- 
tion bands within the present treatment. If we treat the 
correlation effect only between 4f electrons, the number 
of 4f electrons may be rather reduced in comparison with 
the present result. However, the weight of the 4f elec- 
trons is over 2/3 in the total DOS obtained from the ex- 
tracted bands, so that we believe rif vs. U relation may 
not change drastically. 

Figure 4 shows the Sommerfeld coefficient for a few 

7 11 

choices of E { at the very low-temperature region. Note 
that all the obtained coefficients are underestimated 
slightly and systematically since a finite adiabatic con- 
stant 5 is introduced in the numerical calculation (— -> lo + 
id), whose amplitude is given by 5 — 3.0 x 10 -4 Ry. For 
example, 7 is about 12.5 mJ-K _2 -mol _1 for Ej^ 2 = 0.6 
Ry for the non-interacting case, whose value is slightly 
smaller than the result of ref. 2 (13.5 mJ-K~ 2 -mol _1 ). 

With an increase of U, 7 is enhanced and becomes ~ 
30 mJ-K^-mor 1 for U = 0.0175 Ry and E 7 { /2 = 0.6 Ry. 
This result for the Sommerfeld coefficient is roughly con- 
sistent with the experimental one. Taking account of the 
mean valence and the Sommerfeld coefficient observed by 
the experiments, we can estimate the effective Coulomb 
interaction as U ~ 0.017 - 0.019 Ry. 

Let us discuss the temperature dependence of the spe- 
cific heat divided by temperature. We estimate eq. (15) 
by means of the numerical difference, which is defined as 
c v (T). Although the specific heats c v (T) have to vanish 
at T — > for any values of U , the extrapolation values at 
T — > have slightly finite amplitude (Ac v ) due to errors 
on the numerical difference. Thus c v (T)/T is strongly 
enhanced at T —> 0. In order to exclude this artificial 



anomaly, we subtract Ac v from the obtained c v (T) data, 
where Ac v is estimated from c v (T) data at very low- 
temperatures with least square method. The subtracted 
specific heat c v (T) is defined as c v (T) = c v {T) — Ac v . All 
Ac v become very small values, which are of the order of 
1.0 x 10 -2 mJ-K -1 - mol -1 . Hereafter we regard c v (T) as 
the specific heat. 

The figure 5 shows the temperature dependence of 
c v (T)/T, which has to coincide with the Sommerfeld co- 
efficient within the Fermi liquid theory at T = 0. This 
is approximately achieved in our calculation, as is seen 
in Fig. 5. For U — 0, c v (T)/T is almost independent of 
temperature at the low temperature region and corre- 
sponds to the Sommerfeld coefficient. For the interact- 
ing case, c v (T)/T decreases with an increase of temper- 
ature from T = and has the local minimum which is 
smaller than the Sommerfeld coefficient. It seems that 
c v (T)/T increases with an increase of temperature for 
the non-interacting case from T = 0, but decreases for 
the interacting case. With further increase of tempera- 
ture, c v (T)/T has the peak structure in both cases at 
T = Tp Ca k. The schematic temperature dependence is 
shown in Fig. 6. 

Here, we first discuss the behavior of c v (T) /T at T ~ 
for the non-interacting case. It is well known that spe- 
cific heat is proportional to T at low temperature region 
by means of low temperature expansion of the energy 
up to T 2 , so that c v (T)/T is independent of tempera- 
ture. Taking account of the expansion up to T 4 for the 
non-interacting case, the temperature dependences of the 
chemical potential and c v (T)/T are given by 

36 \ D(o) D (o) 2 J ' (ZZ) 
T 3 
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+ 



D (n) = 



TT 4 k 



d n D(w) 



dw r - 



.£,(2) 



D(0) 



T 2 , 



(23) 



(24) 



where Lto stands for the chemical potential at T — and 
D(uj) is the DOS for general systems. 15 For example, 
when the Lorentzian-like DOS D(u) oc (oj 2 + a 2 ) -1 is 
assumed, the second term of the right side in eq.(23) is 
written as 



Il)(2) 

5 



oc - 7a 2 



(25) 



where a (>0) is a half width. When is larger than 
a and positive, is almost filled. Then, this term 

becomes positive. In this situation, c v (T) /T may increase 
monotonically with an increase of temperature from T = 
0. 

The present DOS roughly has a Lorentzian-like struc- 
ture and the large spectrum due to the 4f J = 7/2 level 
below the Fermi level. Since the DOS has a tail at u> ~ 0, 
the similar picture may be realized. Namely, c v (T)/T is 
enhanced with an increase of temperature. Note that the 
above-mentioned behavior occurs in the low-temperature 
region due to the low temperature expansion. 

For the interacting case, while it is difficult to obtain 
the analytical form of c v (T)/T at T ~ 0, the obtained 
result shows that the amplitude of the local minimum in 
c v (T)/T is smaller than that of each Sommerfeld coeffi- 
cient and Cy(T)/T decreases with an increase of temper- 
ature from T = 0. 16 In general, with an increase of the 
Coulomb interaction, the width of the peak in c v (T)/T 
at T ~ becomes narrower since the characteristic en- 
ergy scale is reduced. Thus, c v (T)/T becomes sensitive 
to the temperature at the very low temperature region, 
which is observed in various heavy fermion materials. 

At the low-temperature region (T < 25 K), the sup- 
pression of c v (T) /T of YbALj is observed with an increase 
of temperature from T = 0. The coherence temperature 
(Tcoh ~ 35 K) is one of the typically characteristic en- 
ergy scale of YbAl 3 . Thus, the correlation effect is also 




peak 



Fig. 6. Schematic behavior of c v (T)/T as a function of temper- 
ature. The solid (dashed) line is the expected c v (T)/T for the 
non-interacting (interacting) case. 



reduced by temperature at T <~ T co h, so that the decrease 
of c v (T)/T from T = may occur. 

In the present calculation, although we cannot confirm 
the temperature dependence of c v (T)/T accurately at 
T — > due to the numerical difficulties, we show the 
following possibility; c v (T)/T increases (decreases) from 
T = for the non-interacting (interacting) case at the 
very low temperature region. 

On the other hand, the temperature dependence of 
c v (T)/T has common property for any values of U, in 
which the peak structures appear at T = T pcak . With 
further increase of temperature, c v (T)/T turns to de- 
crease beyond T pea k e.g., T pca k <~ 0.0025 Ry for U = 
and T peak - 0.0007 Ry for U = 0.015 Ry, respectively. 

In the present case, the large DOS exists below the 
Fermi level. Here, the center of gravity of the peak po- 
sition in the DOS is defined as w ~ E pcak . When the 
energy scale of the temperature greatly exceeds the dif- 
ference between -E pea k and Ep, c v (T)/T decays monoton- 
ically as temperature. It is well known that Fermi distri- 
bution function is broadened due to temperature and the 
width roughly corresponds to 4T. The change of gradi- 
ent of the specific heat c v (T) as a function of T occurs at 
T <~ (Ep — -E pea k) / 4, which is shown in Fig. 2 and the in- 
set of 5. For example, each center of gravity of the peak in 
the DOS is roughly located at E peak = -0.015 ~ -0.01 

Ry for U = and at E pcak = -0.006 0.004 Ry for 

U = 0.015 Ry (shown in Fig. 2). Then the gradients of 
the specific heats are changed at T = 0.0025 ~ 0.0035 
Ry and at T ~ 0.0008 Ry, respectively. Furthermore the 
peak temperature of c v (T)/T is rather shifted to lower 
temperature region because of the effect of factor 1/T in 
comparison with the result of specific heat. 

We stress that this peak structure of c v (T)/T results 
from the band structure (or the DOS structure) and the 
peak position is located at less than 1/4 of the difference 
between -E pea k and E-p. With an increase of the interac- 
tion, the effective 4f J = 7/2 level is shifted to the Fermi 
level, so that the peak position of c v (T)/T is also shifted 
to lower temperature region. 

The peak temperature of c v (T)/T is observed at 70 K, 
the amplitude is 80 mJ-K _2 -mol _1 and the Sommerfeld 
coefficient is 40 mJ-K _2 -mol _1 in the experimental re- 
sult. In the present calculation, for U ~ 0.015 — 0.018 Ry, 
the peak temperature of c v (T)/T is located at T pca k ~ 
0.0002 - 0.0006 Ry, which roughly corresponds to 30 - 
90 K. The amplitude of c v (T)/T is range from 40 to 
60 mJ-K _2 -mol _1 . For the value of U to obtain this re- 
sult of c v (T)/T, the Sommerfeld coefficient becomes 20 
- 30 mJ-K- 2 -mol _1 at the very low-temperature region 
(shown in Figs. 4 and 5). Although there are the minor 
disagreements, these obtained results can reproduce the 
experimental ones for YbAl 3 semi-qualitatively. 

4. Summary and discussions 

We have constructed the low-energy effective Hamil- 
tonian with the realistic band structure and discussed 
quantitatively the specific heat for typical heavy fermion 
compound YbAbj. 

The band term consists of the conduction electrons 
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with NEF method, the localized 4f electrons and the 
hybridization term between the conduction and the 4f 
electrons which still keeps the angular dependence of the 
wave-vector k. This model includes only a few parame- 
ters. 

After diagonalizing the band term of the Hamiltonian, 
we extract 14 bands near the Fermi level, which repro- 
duces the LDA result. We further construct the low- 
energy effective Hamiltonian in order to consider the cor- 
relation effect which is mainly caused by the 4f J = 7/2 
electrons. By means of the self-consistent second-order 
perturbation theory with respect to the Coulomb inter- 
action combined with local approximation, the DOS, the 
Sommerfeld coefficient and the specific heat divided by 
temperature are calculated. 

The Sommerfeld coefficient and liniT^o c v (T) /T arc 
also enhanced with an increase of the effective Coulomb 
interaction. The ratio c v (T)/T decreases with an increase 
of temperature for U > from T = 0, but the peak 
structures appear at higher temperature for any U. 

Thus, we stress that the characteristic structures of 
c v (T)/T in YbAi3 may originate from the correlation 
effect and the structure of the non-interacting density 
of states, respectively. Although some disagreements are 
left, the results for U ~ 0.015 ~ 0.018 Ry arc in agree- 
ment with the experiments. 

The peak structure is experimentally observed at ftw ~ 
0.25 eV in the optical conductivity, whose main contri- 
bution comes from the effective 4f J = 7/2 level. While 
the peak structure of c v (T)/T mainly is also caused by 
the effective 4f J = 7/2 level, the corresponding peak po- 
sition becomes rather smaller than the energy difference 
between the effective 4f J — 7/2 level and the Fermi level. 
We can roughly understand the specific heat and the op- 
tical conductivity on equal footing within the present 
effective Hamiltonian. 

In conclusion, the good agreement with the experi- 
mental results reflects the appropriate description of the 
band structure in our model. Thus, our effective Hamil- 
tonian succeeds in describing the realistic heavy-fermion 
materials with conduction bands mainly composed of s- 
and p- electrons. 



Acknowledgement 

The authors would like to thank H. Harima and H. 
Okamura for valuable discussions. This work was sup- 
ported by Grant-in-Aid for Young Scientists (B) (No. 
19740195) from the Ministry of Education, Culture, 
Sports, Science and Technology of Japan. 

1) A. C. Hewson: The Kondo Problem to Heavy Fermions (Uni- 
versity Press, Cambridge, 1993). 

2) H. Kuroiwa, Y. Imai and T. Saso: J. Phys. Soc. Jpn. 76 (2007) 
124704. 

3) H. Harima: Private comunications. 

4) H. Okamura, T. Watanabe, M. Matsunami, T. Nishihara, N. 
Tsujii, T. Ebihara, H. Sugawara, H. Sato, Y. Onuki, Y. Isikawa, 
T. Takabatake and T. Nanba: J. Phys. Soc. Jpn. 76 (2007) 
023703. 

5) T. Ebihara, Y. Inada, M. Murakawa, S. Uji, C. Terakura, T. 
Terashima, E. Yamamoto, Y. Haga, Y. Onuki and H. Harima: 
J. Phys. Soc. Jpn. 69 (2000) 895. 

6) A. L. Cornelius, J. M. Lawrence, T. Ebihara, P. S. Riscborough, 

C. H. Booth, M. F. Hundley, P. G. Pagliuso, J. L. Sarrao, J. D. 
Thompson, M. H. Jung, A. H. Lacerda and G. H. Kwei: Phys. 
Rev. Lett. 88 (2002) 117201. 

7) E. D. Bauer, C. H. Booth, J. M. Lawrence, M. F. Hundley, J. 
L. Sarrao, J. D. Thompson, P. S. Riseborough and T. Ebihara: 
Phys. Rev. B 69 (2004) 125102. 

8) L. H. Tjeng, S.-J. Oh, E.-J. Cho, H.-J. Lin, C. T. Chen, G.-H. 
Gwcon, J.-H. Park, J. W. Allen, T. Suzuki, M. S. Makivic and 

D. L. Cox: Phys. Rev. Lett. 71 (1993) 1419. 

9) S. Suga, A. Sekiyama, S. Imada, A. Shigemoto, A. Yamasaki, 
M. Tsunekawa, C. Dallera, L. Braicovich, T.-L. Lee, O. Sakai, 
T. Ebihara and Y. Onuki: J. Phys. Soc. Jpn. 74 (2005) 2880. 

10) U. Walter, E. Holland-Moritz and Z. Fisk: Phys. Rev. B 43 
(1991) 320. 

11) T. Ebihara, E. D. Bauer, A. L. Cornelius, J. M. Lawrence, N. 
Harrison, J. D. Thompson, J. L. Sarrao, M. F. Hundley and 
S. Uji: Phys. Rev. Lett. 90 (2003) 166404. 

12) H. Okamura, T. Michizawa, T. Nanba and T. Ebihara: J. Phys. 
Soc. Jpn. 73 (2004) 2045. 

13) W. A. Harrison: Electronic Structure and the Properties of 
Solids (Freeman, San Francisco, 1980). Note that the values 
for ro and r c should be interchanged in the solid-state table of 
this book. (See W. A. Harrison, Phys. Rev. B 28 (1983) 550. 
) 

14) W. A. Harrison: Elementary Electronic Structure (World Sci- 
entific, Singapore, 1999). 

15) M. Shimizu: Rep. Prog. Phys. 44 (1981) 329. 

16) Y. Imai and T. Saso: J. Phys.: Conf. Ser. 150 (2009) 042065. 



